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^^ Abstract 
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In this work a study of the gravity is made using Einstein's equation in 
^ , the post-Newtonian approach. This is a method to hnearise the General 

Relativity indicated to treat non-relativistic objects. It enables us to 

construct, from metric-independent elements, fields that are governed by 

^^3 equations similar to the Maxwell's ones in Lorentz gauge. We promediate 

Cn these equations for matter distributed in local systems, like solar systems 

^^ or galaxies. Finally we define the large scale fields for this distribution, 

fj which includes terms analogous to electromagnetic case, like polarization, 

^^ magnetization and superiors terms. 
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1 Introduction 

Due to the fact that are known just two long-range interactions, the gravitation 
and electromagnetism, came the idea there could be some relationship between 
them. The discovery of this link could help for a better understanding of grav- 
ity, especially in the sense of this quantization. Michael Faraday in half of 
nineteenth century, was, apparently, the first to imagine that gravitation and 
electromagnetism have some relationship to each other. Despite his experiments 
have been a fundamental importance for the development of electromagnetism, 
he did not get find any relationship with gravity. Even the negative results of 
the experiments in this direction the Faraday's deep conviction about the exis- 
tence of this connection was not shaken. When Maxwell published the work in 
which showed that electric and magnetic phenomena were manifestations of the 
same field, the electromagnetic, arose the ambition to unify gravity with these 
phenomena, until then described by Newtonian theory. After the reformulation 
of gravity by Einstein the attempts to unify the long-range interactions have 
taken a new direction. In 1918, Weyl [1, published a paper which proposed 
to modify the connection of general relativity by another in which the electro- 
magnetic field would be contained, in the mathematical point of view meant 
leaving the Riemannian geometry. Another attempt to modify the geometry in 
order to unify these two interactions was proposed by Klein and Kaluza, and 
consisted of maintain the Riemannian geometry, but include an extra spatial 
dimension [2Jf3]. The degrees of freedom of gravity in more dimensions, after 
a dimensional reduction, would fall in the usual gravity and electromagnetism. 
Einstein himself was convinced that this relationship should exist and devoted 
the last years of his life to the search of the unified field theory, that is still the 
subject of countless works |:4"}9 . Despite all failed attempts of unification made 
by great names of science persists the belief that this connection exists. 

Historically the attempt to establish an analogy between electromagnetism 



and gravitation began in 1893, when Heaviside 10 investigated how the energy 
propagates in a gravitational field. He proposed a gravitational Poynting vector 
which contained magnetic components of gravitation. In Newton's theory of 
gravitation a problem is the absence of freedom of the gravitational field, which 
is fully defined by a single scalar field, while the electromagnetism by a scalar 
and a vector field. We can associate the Newtonian gravitational potential to 
electromagnetic scalar potential, but to proceed with an analogy between both 
theories we need a vector potential for gravity. In Einstein's theory, the grav- 
itation gains more freedom because now is determined by the 10 independent 
components of the metric tensor. In General Relativity is the gravity who has a 
lot of freedom in relation to electromagnetism. Another problem raised by the 
theory of General Relativity is that its equations are nonlinear while Maxwell's 
are. So, to create an analogy between them it is necessary to linearise General 
Relativity and reduce the degrees of freedom of the metric. Forward was the 
first to write the Einstein's equations in a linearised Maxwell-type structure. 



and proposed experiments to detect the gravitomagnetic fields 11 12 . Bragin- 



ski 13 reached similar equations using the parametrized post-Newtonian ap- 



proximation 14 , while Campbell left in a tensorial form 15 16 . Experiments 



in superconductors was recently proposed by Li 17 19 following the correc 



tions introduced by Ross 20 and deWitt [21] in London equation to include 
the gravitomagnetic fields. 

In this paper we use the post-Newtonian approximation one order beyond 
the Newtonian limit in order to define the gravitoelectromagnetic field with 
independent terms of the metric. The found equations are used to calculate the 
effective fields in large scale, considering the matter distributed in local systems. 
This systems can be planets orbiting a star or galaxies, depending on the desired 
scale. The post-Newtonian approximation is satisfactory in these cases since 



the speed of revolution are much smaller than the light 22 ■ 27 . The deviation 



of effective fields is made in an analogous way which was made by Russakoff 



in electromagnetic case 28 . In this study we expand the contribution of the 
bodies bound to local systems in power series in relation to the size of these 
systems. Finally, we collect the terms analogues to polarization, magnetization 
and superiors terms, such as quadrupole, to define the effective fields in large 
scale. 

This paper is organized as follows: In Sec. II we make a review about the 
post-Newtonian approximation, in particular the expansion of the metric and 
the approximate field equations. We also expand the energy-momentum tensor 
used in this work and obtain general solutions. In Sec. Ill we write the equations 
in the form of Maxwell's equations of electromagnetism, as well as we define the 
gravitoelectromagnetic field with degrees of freedom of the metric. In Sec. IV 
we apply the Maxwell-Einstein's equations for local systems and calculate the 
promediation of the sources making a multipole expansion in terms of typical 
lengths of these systems. In Sec. V we define the effective fields in large scale in 
an analogous way to electromagnetism and write the equations obeyed by these 
fields. In the last section we discuss the conclusion of this work and present 
possibles perspectives. 

2 The post-Newtonian approximation 

To linearise the Einstein's equation we use the post-Newtonian approximation, 
which is a method suitable for objects moving at low speeds compared to the 
speed of light. This approach is recommended to the problem we want to 
study, which is the bodies bounded in solar-type systems. This is because post- 
Newtonian approximation inserts corrections in Newtonian mechanics, and this 
already sufficiently precise that even situations for sending space probes. This 
approach consists in expand the components of the metric tensor and of the 



energy-momentum tensor in powers of the average speed 29 , i.e., we expand 
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where the symbol g denotes the terms of g^i, in order N of typical speed v. 
In this paper we consider one term higher than the Newtonian theory, i.e., we 
will consider up to fourth order in goo, third order in giQ and second order in 
gij . Using the above expansion of the metric in Einstein's field equation we get 
the set of equations 
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where T ^'^ denotes the term of T^" in order of "friv^ /f^, where fh and r are 
the typical values of mass and distance. The above equations been simplified 
by a convenient choice of coordinate system, which obeys the Einstein's gauge, 
defined by 

r^ EE gf^X^^ = 0. (8) 

Since we are interested in this work to apply this formalism to the universe in 
large scale, we can approximate the celestial bodies by point objects. In this 
context we use the energy-momentum tensor of particles that interact gravita- 
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Using the expansion of the metric we obtain the terms needed for the present 
work 
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are the density of matter, the current of matter and the density of kinetic energy, 
respectively. With the components of the energy-momentum tensor we can solve 
the approximated field equations. Considering that the fields vanish at infinity 
and using the Green's theorem [30| we can write the solutions in the form 
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with pg been the gravitational energy density 

Ps = X! '^n(t'nS^{x - X^n). (24) 
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The solution to ^qq was already expected since (f>n is the Newtonian gravitational 
potential, exactly the non-relativistic limit of General Relativity. 

3 The Maxwell-Einstein's equations and the grav- 
itoelectromagnetic field 

In the previous section we find general solutions to the approximate equations. 
In this section we write these equations in a form analogous to Maxwell's equa- 



tions of electromagnetism. Neglecting terms of orders higher than those consid- 
ered in this work, we can write the metric in the form 

ffoo = -1-2(0 + 02), (25) 

g,j = [1-20] 5,,, (26) 

<?« = 4A,, (27) 

where we defined the fields 

= 0„ + V', (28) 

A^ = -X^■ (29) 

In this approach the metric is uniquely determined by four degrees of freedom, 
a scalar field and a vector field A, which satisfy the following equations 

^''^ = ^+4^G[po + P<, + 3pfc], (30) 

V^A = 47rG/o. (31) 

In this approach the gravity has the same number of degrees of freedom that the 
electromagnetic field, so that we can identify them as four components of the 
electromagnetic quadrivector A^ . Another indication that we are in a correct 
way to build an analogy between gravitation and electromagnetism is that the 
not trivially satisfied component of harmonic coordinates fixation ([8|, gives, in 
this approximation, the condition 

^+V.A = 0. (32) 

That, analogous to electromagnetism, is the Lorentz gauge. This condition 
allows us to define the gravitoelectric and gravitomagnetic fields analogously to 
electromagnetism, i.e. 

3 A 

Bg = V X A (34) 

With the same calculation used in electromagnetism we obtain, in this approx- 
imation, the Maxwell-Einstein's equations 

ipk] , (35) 

(36) 
(37) 

(38) 



V.Eg = 


= -—[pa + Pg 

-- 0, 


V y^Eg = 


dBg 
dt ' 


yy.Bg = 


-, , dEg 
-- P9JO+ ., 



where we define the constants 
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Analogously we can obtain the equations for the fields separately, which are 
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In the last equation was neglected the term with second derivatives of Bg as 
that term is higher order them the approximation used here. We can use the 



Green's theorem to obtain the solution of (40) and (41) which are 
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where the bracket [ ]ret rneans that t' is the retarded time, t' — t — \x — x'\. 

The differences between these equations to Maxwell's is the negative sig- 
nals in sources of gravitoelectric, and gravitomagnetic, fields which gives due to 
the fact that gravity is attractive while electromagnetism is repulsive for equal 
charges. Another difference is that the source of gravitoelectric field is not only 
the matter density but contains terms proportional to the density of gravita- 
tional energy and kinetic energy, which comes from the fact that not only the 
mass interacts gravitationally but also the energy. Although we had show be 
possible to construct fields analogous to electromagnetic with the independent 
terms of the metric the analogy can not be taken later on. This is because is 
not possible write the equation of motion in the form of the Lorentz force. The 
geodesic equation, in this approximation may be written as 
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which has terms in direction of velocity, and bilinear terms, which does not 
allow us to write it in the form of Lorentz force. 



4 Deviation of large scale equations 

In this section we derive the effective Maxwell-Einstein's equations in large-scale, 
analogously to the electromagnetic case 28p0 . In electromagnetic case is made 




Figure 1: Local coordinate systems. It is used to write the position of the bodies 
connected to local systems, where a;„ orients the center of mass of n-th local 
system and Xj the j-th bounded body, respecting the relation Xj — Xn + Xjn- 



a power series expansion on the molecules internal distances, considered much 
smaller than the scale of interest. In this work we will do a similar procedure, 
but we will expand on the internal distances of local systems. These systems 
can be solar systems or galaxies, since they are similar bounded structures, 
depending on the distances of interest. For this we will use a coordinate system 
similar to that used by Russakoff 28 , which localize the bodies bounded in 



local systems in relation to the mass center of the system, as illustrated in FIG. 
[T] Considering that the bodies involved produces fields that obey the equations 
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we must compute the spatial average of the amounts above to obtain the effective 
equations for larger scales. The spatial average of a function F{x,t) in relation 
to a function test w{x) is defined as 



{F{x,t)) = / d\'wix)F{x~x,t), 



(49) 



where w{x) is a real function, non-zero in a neighborhood of :? = and normal- 



ized to unity in the whole space 28 . The test function need not be specified in 
detail, all that is needed are general properties of continuity and smooth vari- 
ation. This allows the rapid convergence of the Taylor series development over 
distances of the local systems size, i.e., that can write 

w{x — Xn — Xjn) ~ w{x — Xn) ~ Xjn-Vw{x — Xn) + 

+ iv.V.[f,„f,„t.(f-f„)]. (50) 

Respected the condition for the test function we can define the large-scale 
gravitoelectric and gravitomagnetic fields. Eg and Bg, as the averages of fields 
Bg and bg, i.e., 

Eg =<eg> e Bg =< bg >, (Sl) 

so, the promediate Maxwell-Einstein's equations are 

^■Eg - -^[{po)+3{pk) + {pg)], (52) 

V.Bg = 0, (53) 

Vx4 = -^, (54) 

yy<Bg = -pgQo) + ^- (55) 

To find the desired equations we need to calculate the spatial averages of po, Pg, 
Pk and Jo J that we will calculate below. 

Calculation of spatial average of pQ 

Considering that in large-scale the involved bodies can be taken as point, in the 
coordinate system shown in FIG.jT] the density of matter, po, can be written as 

Po = X! '^J^i^ ~ ^j) + X! X! "^j^i^ - Xn- Xjn), (56) 

j(free) n j(n) 

where Xn orients the center of mass of the n-th local system. The first term 
of Po is the contribution of free bodies while the second term represents the 
contribution of these systems. The average of density of matter of n-th local 
system is, using the definition (|49|), given by 



< Po >n=^^rnjw{x- Xn- Xj„). (57) 



As we are considering the dimensions of these local systems much smaller than 



the dimensions considered, we can use the expansion ( 50 ) to write 



+ - ^ mj(f-,-„.V)[fj„.Vw(x - f„)]. (58) 

As Xn guides the center of mass of the n-th local system, then J2j(n) ^^j^jn — 0. 
So we can write 

< Po >«~< mnd{x~ Xn) > +V. [V. < q„(5(f - f„) >] , (59) 

where 

q« = 2 51 ^j^jnXjn, (60) 

Jin) 



is the quadrupole moment of the n-th local system. Finally, according to (56), 
we can write the average of density of matter of the form 

<Po>«P/ + V.[V.Q], (61) 

where we define the density of free matter and the quadrupole momentum 

Pf = ( ^ mj5{x ~ Xj)) + (^ mnS{x ~ f„) ) , (62) 

\jifree) / \ " / 

Q=(^q„<5(f~x„)). (63) 



The result (61) is similar to obtained in electromagnetic case with electric charge 
exchanged by mass. This result has not analogous term to dipole momentum 
as that term represents the sum of the momentum of each local system, which 
vanish due to reference system taken on center of mass. 



Calculation of spatial average of jo 

Analogously to what was done in calculation of average of po we can write, using 
the same coordinate system, jo as follows 

jo = Yl mjVj5{x-Xj) + 

Jifree) 

-|-^^mj(z/„ +Vjn)S{x- Xn - Xjn), (64) 

n j{n) 



dx^ 



where Vj — -^ . Again, the first term is the contribution of free mass while the 
second is the contribution of the bounded bodies. We will calculate separately 
the average of bounded term which we can expand in series. The average of 
current matter of the n-th local system is given by 



< jo >ri = ( ^ mj {Vn + Vjn)5{x - X„ - Xjn) ) 
\j(n) / 

W mnVnW{x - Xn) + ^n V. V. [q„ w(x - X„)] - 

-^{mjVjnXjn)-'^w{x- Xn), (65) 

where was neglected terms proportional to VjnXjnXjn- To simplify the above 
equation we use the following identity 

^(mjWj„fj„).Vw(f - x„) = — [q„.Vw(f - f„)] + 
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so that we can write the average of current of matter of n-th local system ( 65 1 
in the form 

d 

< Jo >n ~ {m,iVnS{x ~ Xn)) ~ ^V. (q„w(f - f„)) - 
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+V X {jln5{x-Xn)) , 



where 



l^n — ^ / ^ fTT'j [Xjn X Vjnjj 
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is the gravitomagnetic momentum of n-th local system. Finally, according to 
(64 ), we can write 



Jo 



where we defined 
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(69) 



jf = ( X! "f^J^'J^i.^ ~ ^i) ) + ( X! '^n^^"^(:? - Xn) ) , (70) 

\jifree) / \ n / 
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M= (^„(5(f-x„)), (71) 

as the free current density and gravitomagnetization. The result is similar to 
that achieved in the electromagnetic case, with exception the polarization term, 
which vanish in the gravitational case. 



Calculation of spatial average of pk 

We will now calculate the spatial average of the term proportional to pk which 
represents the contribution of kinetic energy. Similarly to what was done in the 
previous promediation we will write this term as 

jifree) 

+ ^'^'^mk{Vn+Vkn) S{x-Xn- Xkn)- (72) 

n k(n) 

Again, the first term is the the kinetic energy of free bodies while the second is 
that of objects bounded in local systems. Let us look separately the contribution 
of the bounded term 

{Pkjn = ^\^^'nnk{Vn+ Vknf 5{x - f„ - Xfe„) \ 
\k(n) I 

~ ^rrinvlwix ~ Xn) + -V .V .qnvlwix - Xn) - 



Vn- ^[mfeWfe„ffe„].Vw(x - f„) + 
fc(n) 

- Y^ nikvl^wix - Xn). (73) 



k{n) 

To improve this expression we use the identity 

Vn- ^[mfeWfc„Xfe„].Vw(f - x„) = V.V. [qnvlw{x - Xn)] + 

k(n) 
+ — V. [c\ri-VnW{x - Xn)] + V. [Vn X ^„w(x- f„)] - 

at 



-V. 



q„.— w(x-x„) 



+ V.V.[(q„ X ■?„) X u„w(f-x„)] , (74) 
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to write the eq. ( 73 ) as 



(Pfe)n ~ -(m„i;^(5(f-f„)) - -V.V.(q„w,2,(5(f-f„)) + 



9 ( Yl ™kvlJ{3 



2 , 

\k(n) I 

u / dv 

- — V. (q„.v„(5(f - f„)) + V. ( q„.-^(5(f- f„) 

-V.V. ((q„ X w„) X w„(5(f - f„)) - V. (u„ x ^„(5(f - x„)) .(75) 

As each body follows the Newtonian equation of motion |29j , we can write the 
following term as 



m;ra^n 



^m ^n| 



The first term represents the contribution of free bodies and can usually be 
ignored, but will be kept here for completeness, while the second term represents 
the contribution of local systems, which had been expanded in power series. 



Finally we can return to equation ( 72 ) to obtain 

(Pfe)«(p/c)/-V.P,.-V.V.Qfc, (77) 

where the free kinetic energy was defined as 



\]{free) / \ n / 

+ o{YY^kvlj{x-Xn)) , (78) 



\ n k[n) I 

the kinetic polarization and the kinetic quadrupole momentum as 

Pk = {y^Vn y- fln^ix - Xn)) - (^ q.n-Gmm—^ ^(5(x-f„) 



m •^n I 



Ju 'I Ju J- 



-(J2 Yl q»-GTOj- ^^ " S{x-Xn)) + 
\ n Jifree) I^J ^"1 / 

Qfe = 9 ( X!^""""^*^^^ ^") ) + ( X!(^" ^ ^") ^ w„(5(f - f„) \ . (80) 
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This term has no electromagnetic equivalent but behaves similarly to the average 
of density of matter. Unlike previous promediations that presents polarization 
terms. 



Calculation of spatial average of pg 

Finally we calculate the average of term that represents the contribution of 
the Newtonian gravitational energy in source of gravitoelectric field. Using the 
coordinate system defined in the average of density of matter, FIG. [T] we can 
write it as 
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(81) 



The first term above is the gravitational interaction of free bodies. The second 
and third is the interaction between free bodies and bounded bodies, while the 
last is the interaction between bounded bodies. Again, the terms summed on 
the free bodies can be neglected, but we will keep these terms for completeness. 
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Let us look separately the average of the last term 

\ n j.k ■' I 

+V. (2 2^ q™. — —3 d{x - Xrn) ) - 

\ m.n F™ ~ -^"1 / 

-V. ( 2^ 2^ I ^ _ 1 ' (5(x -Xn)) + 

\ n j,k I^J" ^*="l / 

\ n j,k ' -^ ' I 

+^-^-(e|^'^(- --»))• (82) 

where in the last passage we expand \xm — x„ — (xkn — Xjm)\ around Xkn — 
Xjm = 0, for m ^ n, w(x — Xm — Xjm) around Xjm = and we eliminated higher 
order terms. Now we focus on the average of the third term of pg, that is 

\k(free) rn j[rn) ' "^ ' / 

\k(free) m j(m) I J *= I / 

-v,(i:i:i:^^r=fp.-..). 

\k(free) m j(rn) l"^-? "^'=1 / 

+V.V./- y y- y- ^'".^"^^•^i-f'^'" J(f - a-^)\ ■ (83) 

Finally, adding the contribution of all terms we can write the average of pg as 
< 99 >~ (Ps)/ + V.Pg + V.V.Q^, (84) 
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where we define the free gravitational energy 

\j(/ree) " fc(„) I J *= I / 

\fc(/ree) "i j(m) I J *= I / 

\ . , \^ Tn '^kn 1 

\ n j^k •' I 

and the quantity 

5 o / Y^ Gm„(a?m-a?„) ^ \ 

"g = 2( 2^q„i. — — —3 — d(x - Xm) ) - 

\ m n F"i ^ 2;„| / 

~ \ 2^ 2^ I ^ _t , six - x„) ) - 
\ « j.fe i^j" ^'="1 / 

\k(free) ra j(m) I^J ^'^l / 

Qs = (22^ 2^2^ — \x--k\ '^(^'^-)) + 

\ k(free) ni j{m) \ J f^\ / 



m,n 



2^ 1^ _ ^ I "^ "'' / ^ 



are the gravitational polarization and the quadrupole monicntum of gravity. 
Again, despite having no electromagnetic corresponding, the average this term 
exhibits behaviour similar to other sources of gravitoelectric field, including 
polarization terms. 
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5 Definition of Large Scale Fields 

In the previous sections we calculated the averages needed to write the prome- 
diated Maxwell-Einstein's equation. Finally we find the result 

Pf+'^{Pk)f + iPg)f + 



V.S,«- 



P3-3Pfc + V.(Q-3Qfc + Qg) 



V.S„ = 0, 



V X Ea 



V X Ba 



dt ' 

-PgJf+Pg-Ql'^-Q 



-fIgV X 



M - V. / ^ q« X Vn6{a 



dt 



The equation ( 88 ) leads us to define the field and the effective density 

1 



D„ 



En 



Fg-3Pfc + V.(Q-3Qfc + Qg) 



Pe = Pf +i{pk)f + {Pg)f, 

in order to validate the equation 

V.D„ 



1 



-Pe 



The definition (92) allows us to write eq. (91 ) in the form 

dD 

where we define the effective current and the large scale field 
^ - d 



Je 



3f 



dt 



Pg-3Pfc+V.(Qg-3Qfc) 



Hg^ Bg+ Pg\/ X 



M-v.{J2 q™ ^ ^nSi 






(88) 
(89) 

(90) 

(91) 

(92) 
(93) 

(94) 

(95) 

(96) 

(97) 



According to the above definitions, we can write the large scale Maxwell-Einstein's 
equation in the form 

1 



V.i^„ 



V.P„ = 0, 



V X Ea 



dBg 

dt 



\7 X Hg = -Pgje 



dDg 

dt 



(98) 
(99) 

(100) 

(101) 
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Finally, we can note that both the free and the effective mass are conserved, i.e. 

= V.// + ^, (102) 

= V./e + ^. (103) 



The conservation law ( 103 1 indicates that the effective current is a current as- 



sociated with source of gravitoelectric field. This indicates a symmetry in large 



scale equations much greater than in fundamental equations (35) - (38). In 
latter the source of gravitomagnetic field was only the flow of matter, without 
regard the flow of other forms of energy, whereas the source of the gravitoelec- 
tric field take them into account. In large scale equations that asymmetry was 
broken, because the source of the field Hg is the same order flow of the source 
of field Dg. 

6 Conclusions and Perspectives 

We show in this paper that it is possible to construct fields with independent 
elements of the metric in the fourth order post-Newtonian approximation, which 
obey equations similar to Maxwell's electromagnetism. Through these equations 
becomes explicit that not only the matter interacts gravitationally also more 
energy - in this case the gravitational potential and kinetic energy. Due to gauge 
obeyed by these fields, from fixing coordinates, the propagation of these fields 
occurs at the speed of light, which agrees with the spirit of special relativity. 
Despite this progress in building an analogy between both theories, it is unable 
to write the equation of motion in the form of the Lorentz equation. With 
gravitoelectromagnetic fields we found the effective fields for a distribution of 
matter organized in local systems, as occurs in the universe where matter is 
organized into bounded structures as solar systems, on a scale, galaxies, on 
a scale greater. Only after this last level is that we can use a cosmological 
model, as Friedmann, which presupposes matter homogeneously distributed. 
The equations for these effective fields, Dg and Hg, exhibited a greater symmetry 
than the equations of fundamental fields, since their sources satisfy an equation 
of continuity in the same order, i.e., the source of field Hg is the flow associated 
with the source of field Dg. We show that the definition of these fields are 
greatly simplified by choosing the coordinate system fixed at the center of mass 
of local systems, where the terms of dipole vanish. As perspectives to further 
studies to be made based on the results obtained in this work we can mention 
the inclusion of the cosmological constant, as well as calculate the corrections 
of order higher than those considered in this work. 
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